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Abstract 

We extend hadronic models for p-meson propagation in cold nuclear matter via 
coupling to in-medium pions to include finite three-momentum. Special care is taken 
to preserve gauge invariance. Consequences for photoabsorption on the proton and 
on nuclei as well as for the dilepton production in relativistic heavy-ion collisions 
are discussed. 



1 Introduction 



The investigation of matter under extreme conditions and the modifications of hadron 
properties with density or temperature is one of the main topics in intermediate and 
high-energy nuclear physics. Experimentally, the cleanest information is obtained from 
electromagnetic probes which penetrate the medium almost undisturbed. In the vector 
dominance model (VDM) Q electromagnetic processes are mediated by neutral vector 
mesons [p^, u and </>). Therefore the medium modifications of vector mesons are of par- 
ticular interest. For the p meson there are both, theoretical and experimental indications 
that its properties are changed in nuclear matter: 

Using scale invariance arguments. Brown and Rho concluded that the mass of the 
p meson and other hadron masses should drop as a function of temperature or density. In 
fact, many authors predicted dropping vector meson masses, e.g. Saito et al. within the 
Walecka model or the Guichon model . From a QCD sum rule analysis Hatsuda and 
Lee also concluded that the mass of the p meson is shifted downwards considerably 
at nuclear matter density. On the other hand, as recently shown by Leupold et al. 0, 
predictions from QCD sum rules are not unique unless further assumptions concerning the 
p meson width are made. Klingl et al. ^ could satisfy the QCD sum rules with a model 
which predicts an increased width of the p meson while its mass stays almost unchanged. 

A similar picture emerges from the dilepton data measured by the CERES collabora- 
tion in ultra-relativistic nucleus-nucleus collisions : the strong enhancement of dilepton 
production at low invariant masses (~ 0.3 - 0.6 GeV) was interpreted by Li et al. ||T^ as 
a signature for a dropping p-meson mass, supporting the scaling hypothesis of Brown and 
Rho 1^. Alternatively, however, the CERES dilepton spectra could also be explained by 



a strong broadening of the p meson in hot and dense matter [11, O 



There are two classes of "conventional" processes which potentially cause a broadening 



of the p meson: First, as was demonstrated in the models of Chanfray and Schuck [|I3 



and of Herrmann et al. [|14[ for cold nuclear matter, the pion cloud, which gives rise to 



the vacuum width of the p meson, is modified through interactions with the surround- 
ing medium. A second important effect is expected from baryonic resonance formation 
through direct interactions of the p meson with nucleons. This was considered first by 
Friman and Pirner and investigated extensively by Peters et al. ||16|| . 



Combining both mechanisms, extended to finite temperatures, Rapp et al. fl^ ob- 
tained a reasonable description of the CERES data. However, the underlying model for 
the pionic part [0 was valid only for p mesons at rest {q = 0). Therefore it was as- 
sumed in ref. that this part of the selfenergy does only depend on the invariant mass 
M of the p meson, but not on its 3-momentum q. For the resonance contributions the 
g-dependence was taken into account. 

The aim of the present paper is to construct a model for the 3-momentum dependence 
of the pionic part of the p-meson selfenergy at zero temperature. To a large extend we will 
follow the work of Herrmann et al. [14]. Our article is organized as follows: in section ^ 
we give a short summary of the VDM description of the p meson in vacuum. In section ^ 
the modification of the pion propagator in dense matter is discussed. The corresponding 
corrections to the pmr and ppTTir vertices, as required by Ward-Takahashi identities, will 
be derived in section Using these vertex functions and the in-medium pion propagator 
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we construct the p-meson propagator in matter as detailed in section ^. In section ^ we 
discuss two applications of the model. First we calculate photoabsorption cross sections 
on nucleons and nuclei and show how these processes can be used to constrain some of the 
model parameters, like the tcNN form factor. Since real photons have a fixed invariant 
mass M = 0, the 3-momentum dependence is necessary to obtain nonzero contributions 
to the cross section. Finally we will discuss the impact on dilepton production rates in 
heavy-ion collisions. 



2 The p Meson in Vacuum 

The p meson has a hadronic width Fp 150 MeV with two-pion decay accounting for 
~ 100% of it. We restrict ourselves to the neutral p meson, i.e. tt+tt" decay. The free 



Lagrangian involves the isovector pion field and the neutral vector field p^ and reads [|17 



C^ + C, = ■ d^$~ ^ml$ -^p^^p^^ + l{mffp^p^ , (2.1) 

where p^y = d^p^ — d^p^ is the field strength tensor of the p field, and m!f'> the bare p 
mass. Minimal substitution, 

dJ-^{d^ + igp^T;)$ , (2.2) 
[g = vrp-coupling constant) leads to the interaction Lagrangian 

C^p = \i9P,{T:,4>- 9^0 + S'^^ ■ T30) - i(7VP^^30 ■ T30 , (2.3) 

which contains pixTC and pp-R-K vertices. 

To second order in (7 the p- meson selfenergy in vacuum is given by 



-iSp^(g) = 9^ J 



d^k {2k + q)^{2k + q). 



{2tiY {{k + qf -ml + ie){k^ - ml + ie) 

The corresponding diagrams are shown in fig. 1. 

In the VDM the p meson couples to a conserved current. Consequently, the selfenergy 
must vanish for = and has to be 4-dimensionally transverse: 

g^Sp, = g'^Sp, = . (2.5) 

This is in general not fulfilled if the divergent integrals in eq. (|2.4|) are regularized by a form 
factor as was done in ref. |13[. Therefore we will employ the Pauli-Villars regularization 
scheme which is known to preserve gauge invariance and hence transversality. As each 
integral diverges quadratically, we add two regulator terms, 

2 

S^;.(g; m^) — ^ ^tMuiq; m^) + ^^^^^(g; M^) . (2.6) 

i=l 
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Figure 1: p-meson self energy in vacuum. 



Here we deviate from ref . [0 , wliere first botli integrals are added, and then the remaining 
logarithmic divergence is canceled by only one regulator. Instead we choose 

ci = -2, Ml = ^ml + A2 , 
C2 = 1, M2 = y< + 2A2, 

(2.7) 

with one free cutoff parameter Ap which is taken as 1 GeV, a typical scale for hadronic 
models. 

Iterating the selfenergy insertions in a Dyson equation by standard techniques the full 
p propagator of the vacuum model is obtained. The two remaining parameters, the vrp 
coupling constant g and the bare p mass m'^p^ are fitted to the pion electromagnetic form 
factor. For g = 5.9 and mf-p^ = 853 MeV we obtain the fit shown in the left panel of fig. 2. 
With these values the p-wave nn scattering phase shifts Si are also reproduced reasonably 
well (right panel of fig. 2). 



3 The Pion in Nuclear Matter 

In nuclear matter the interaction of a pion with surrounding nucleons generates a pion 
selfenergy. The most important contributions are particle-hole [Nh) and Delta-hole {Ah) 
excitations arising from p-wave pion-nucleon interactions. In the present work they will 
be treated in non-relativistic approximation. However, it will turn out to be useful for 
the systematic derivation of the coupling to the p meson (see next section) to start from 
relativistic Lagrangians and then take the non-relativistic limit. 



For the vrA^ interaction, we use 17, 21 : 



Cn = i){i<^ - mN)i) , (3.1) 
C-.N = —tpi'Yr^-dJ, (3.2) 

with a TrA^-coupling constant /^/(47r) = 0.08. The standard non-relativistic Feynman 
rules can be derived formally from these Lagrangians by neglecting antiparticle con- 
tributions and expanding relativistic vertices and spinors to lowest non-vanishing or- 
der in l/rriN. However, the relativistic kinematics in the particle and hole propagators, 

^n{p) = \jv^ + "^AT, will be kept. 
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Figure 2: Left panel: Electromagnetic form factor of the pion as a function of 

in GeV^ . The data are taken from refs. [18, [7^ . Right panel: vrvr scattering phase shift 5j 

: GeF. The data are taken from 



in degrees as a function of the invariant mass M^^ 
ref J§/. 



For a relativistic treatment of the spin-3/2 Delta resonance (A), the formahsm of 
Rarita and Schwinger |^ can be apphed. The free A Lagrangian and the ttA^A interaction 



Lagrangian read [21, 22 



Ca = -i^A"^^ - m^W + + l.d^W - + mA)lur , (3.3) 



NA = - — 0H^,-dy + h.c. 



(3.4) 



For the vrA^A-couphng constant we adopt the Chew-Low value /a = 2/7V. The vertices 
derived from C.„na and the Rarita- Schwinger spinors are expanded in powers of 1/mA, 
and again only the lowest non-vanishing order is retained, leading to standard Feynman 
rules. As for the nucleon propagator, we keep relativistic kinematics, 00 /^{p) = 



p2 + rn\, 

but neglect the antiparticle contribution to the A propagator. To account for the finite A 
width. Fa ~ 120 MeV, we introduce a constant imaginary part in the denominator of the 
A propagator. As will be seen later, inclusion of the momentum and energy dependence 
of the delta width would lead to enormous complications in maintaining gauge invariance 
of the p-meson selfenergy. 

In this approximation the Nh and Ah contributions to the pion selfenergy are given 

by 

{^Nh,Ah)ab{k) = k'^SabTiNh,Ah{k) . (3.5) 

with isospin indices a and b, and the dimensionless Lindhard functions 
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^m^^ J\p\<PF {2iTy 



1 



- ujn{p + k) + ujn{j>) +ie P + oj^ip + k^) - ujn{j)) - ie 



(3.6) 



and 



16//a\2 f d^p 



9 ^m^^ -'|p1<pf (2vr)-^ 

( \ I ) 

^fco - wa(p + k) + ujn{p) + /cO + o;a(p + fc) - cu^(p) - '^vj 

(3.7) 

Because of the relativistic kinematics, only the angular integrations can be evaluated 
analytically, whereas the integration over |p| must be performed numerically. 

Phenomenology furthermore requires to include repulsive short-range correlations in 
the particle-hole bubbles, parametrized in terms of Migdal parameters g'-^^ {NNNN ver- 
tex), g'12 {NNNA vertices) and (?22 {NN AA vertex) |^. These interactions lead to a 
renormalization of the ttNN and ttNA vertices and induce a mixing of the Nh and Ah 
excitations: 

j-j _ IlAffe + IlAfe - {g'li - 2g[2 + g22)^Nh'n^Ah ^-^ 
1 ~ g'u^Nh — 922^Ah + (911922 ~ 9l2)^Nh^Ah 

Finally we introduce a monopole form factor at the tcNN and ttNA vertex which, in our 
non-relativistic approximation, depends only on the three-momentum k of the pion, [| 

r.(^) = . (3.9) 

A2 + k^ 

The final result for the pion selfenergy thus becomes 

S;,(A;) = k'6abll'{k) = kXbTUmik) , (3.10) 
and the in-medium pion propagator is obtained as 

GAk) = J = / r.^,.,. • (3-11) 



mi 



S'(A;) k^-ml- k^'ik) 



To render the numerical (in-medium) calculations in sec. ^ more tractable most of our 
calculations will be performed within the so-called "3-level model" . In this approximation, 
which is equivalent to neglecting the Fermi motion of the nucleons, the pion propagator 



^ Note that this form factor is different from (A^ — toJ)/(A^ + k'^), used in re fs. [fl] , [I^ , which 



was adopted from the Bonn potential [g5|. Alternatively one can choose the form (3^) with a redefined 
coupling constant fN{A){k'^ = 0) |^^. This choice is more appropriate in the present context since 
^Tr{k — 0) is normalized to unity, independent of A. 



6 



is a mixture of three quasiparticle propagators. Besides reducing the numerical effort it 
allows for a more transparent physical interpretation of the results. 



Let us briefly outline the main features of the 3- level model. We start from eqs. (3.6) 
and (|3.7|) for the Nh and A/i polarization. Assuming the momentum |p| of the hole to be 
small compared to the pion momentum |/c|, one obtains 

n^.(fc)= ^ "f^!^ . , n^,{k)= /^^fL , (3.12) 

^0 ~" ^Nh{k) + i£ — ^Ahik) 

where 

^Nh{k) = ujN{k) - niN , aNh{k) = A^i—) ^Nh{k) , 

^Ah{k) = ujA{k) -niN - -Fa , aAfe(^) = — 7r^(—) ^Ah{k) . 

2 9 37? ^m^^ 

(3.13) 

Substituting expressions ( p.l2| ) into eq. ( ^^ , the pion selfenergy can be written as 

m) = + , (3.14) 

with «! 2 and i7^2 being some combination of aAr/i, oa/i, ^at/i and ^^A/i- Again, this result 
must be multiplied by F^: 

nW = ^^^ + ^^^, (3.15) 



where q;']^ 2(^) = r^(A;) cti 2(fc). Finally we insert expression (|3.15|) into eq. ( p.ll|) , and 
obtain a pion propagator of the following structure: 

G,W = ^l<*L + ^iW^ + ^lW^. ,3.16) 

ki- ^Kk) ki - ^i(k) *•§ - uiik) 

The pion is a mixture of three quasiparticles with dispersion relations wi, uj2 and tus, i.e. 
the so-called (A^/i)i-branch, the pion-branch and the ( A/i)L-branch, respectively. Here 
the index L means "longitudinal" , indicating that only spin-longitudinal excitations (i.e. 
(X a ■ k or S ■ k) can mix with the pion. On the other hand, Qi and fl2 can be identified 
as the dispersion relations for spin-transverse excitations {Nh)^ and {Ah)^ which do not 
mix with the pion. Fig. 3 shows these three branches for g = qq = O.IQ fm~^ as well as the 
dispersion relations Qi, uj^^ and VL2 of non-interacting Nh, pion and A/i (quasi-)particles, 
respectively. The corresponding strength functions 5*1, 5*2 and 5*3 are displayed in the 
right panel of fig. 3. 
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Figure 3: Pion properties in the 3-level model calculated at nuclear matter density Q = Qo = 
0.16 fm^^ using A = 1200 MeV, g'^ = 0.8 and g'12 = §22 = 0-5 P^, W^- Left panel: Dispersion 



relations of the various branches. Solid lines: a;i^2,3 (from bottom to top), long-dashed lines: 
0,1^2- The short-dashed line indicates the dispersion relation of the free pion, uj-,^. Right panel: 
Strengths of the three branches. Si (solid line), S2 (long-dashed line) and S3 (short-dashed line). 

4 pnn and ppmr Vertex Corrections 

The transversality of the p-meson selfenergy can be inferred from Ward-Takahashi iden- 



tities which relate the pTrvr and ppmr vertex functions to the inverse pion propagator |14]: 



q'K^!iik,q) = ge,,,(G-\k + q) - G-\k)) , (4.1) 
q^T'^tlik, k, q) = ig{e,,j:'}!}{K -q) - e,,r}!^{k + g, -q)) . (4.2) 

In eq. ( f4.2|) we have restricted ourselves to the case of equal in- and outgoing pion mo- 
menta {ki = /C2) which is sufficient for our purposes. The vertex functions r|^^^^^(fc, g) and 

^'ilulbiki, k2, q) are illustrated in fig. 4. 

It is convenient to split the full vertex functions into a bare part and a vertex correction: 



r;S(A:,g) = ge,ab{2k + q), + f'^2{Kq) , 
^'ilulbiki,k2,q) = 2ig'^{5ab-hahb)gt,v + '^'^tlb{ki,k2,q) . (4.3) 

In the vacuum the Ward-Takahshi identities are trivially satisfied. In nuclear matter, 
where the pion propagator is modified as discussed in the previous section, vertex cor- 
rections must be taken into account. They can be constructed by drawing all possible 
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a 



a 




Figure 4: Vertex functions r'^^l{k,q) (left) and T'^^li^{ki,k2,q) (right). 

selfenergy insertions for the pion propagator and then couple the p meson to all possible 
lines and vertices. 

As a first step we will calculate the vertex corrections F^'^^'^) and r('*^'^) which corre- 
spond to the pion selfenergy T^Nh, i-e. Nh excitations without the ttNN form factor. For 
this purpose we have to specify the interaction of the p meson with nucleons. To ensure 
gauge invariance, the interaction Lagrangians are constructed from the Lagrangians C^, 
jO-ttN (eqs. (p. ID and (|3.2|)) by minimal substitutions similar to eq. (^.21): 



^PN = -|#r3^, (4.4) 

jO-pnN = ig—ip'y^/>Ti>-Ts(f). (4.5) 

For the resulting vertices we again perform a non-relativistic expansion keeping only the 
leading order (l/mAr)°. The relativistic pNN vertex and its non-relativistic expansion 

rr~' = -»V3-l (4.6) 

Thus, to order {l/m]\i)^, the p meson couples only to the nucleon charge and not to the 
convection- or magnetization currents. In fact, one could have added a term to eq. ( [4.4| ) 
which describes the tensor coupling of the p meson to the anomalous magnetic moment 
of the nucleon and which is gauge invariant by itself. However this term would also not 
contribute to order (l/miv)°. 

The relativistic pNN vertex fulfills the Ward-Takahashi identity 

q>^Tl^^^)(p,q) = -^gIl^=-^g^(^S^\p + q)-S^\p)) , (4.7) 

where p is the nucleon momentum and q the momentum of the p meson. In a consistent 
non-relativistic approximation an analogous relation should hold if we replace the rela- 
tivistic nucleon propagator Sp by the non-relativistic one, Gj\f. However, starting from the 
r.h.s. of eq. (f4.6|) the Ward-Takahashi identity is violated at order (l/rriAr)^. To correct 
for this, we simply add the missing term to the 0-component 

" O(^) /i= 1,2,3 . 
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Figure 5: pvrvr vertex correction diagrams due to Nh-loops. 

At least for time-like momenta, q, this approximation is of the same accuracy as the non- 
relativistic expansion. A similar approximation is necessary for the following expression 
appearing in some ppnn vertex corrections, 

(4.9) 

Using relativistic Feynman rules, i.e. Gat Sp and F^''^^^ — > —ig^'^^-, we find 

?i ?2^/.<^ = -^-r (Sf{p) + Sf{p + qi- q2) - Sf{p + qi) - Sf{p - ^2)) • (4.10) 



4 

Within our non- relativistic approximation, the corresponding relation (with G^ instead 
of Sp and F|f^^^ according to eq. ( ^78|) ) is again violated at order (l/mjv)^. We correct 



for this again by simply adding the missing term to the 00-component: 



^ ^) -i^^{GN{p) + GN{p + qi-q2)-GN{p + qi)-GN{p-q2)) fi = u = , 
otherwise . 

(4.11) 

From CpT,N we deduce the pirNN vertex as 

Yip^NN), ^ , Ik 5 , ^ I ° + ^(^) ^ = ° ' (4.12) 

\ -^9i^es,aa^n + 0{^) /i = 1, 2, 3 . 

We are now in place to calculate the vertex corrections corresponding to the Nh 
selfenergy of the pion. Two of the four possible pTcrc vertex correction diagrams are shown 
in fig. 5, the remaining ones are obtained by interchanging the pion lines. For the ppmr 
vertex corrections, there exist sixteen possibilities to attach two p-meson lines to the Nh 
bubble. Eight of these diagrams are shown in fig. 6 and the others can be generated by 
interchanging the pion lines. 

Using eqs. ( |4.8| ) and ( 4.12| ) and the definition ( p.6|) , we obtain the following correction 



for the pvTTr vertex function 



^{3 Nh). 



^iab \k,q) = ge3ab {{k + q)iUNhik + q) + kiIlNh{k) 



(4.13) 
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Figure 6: ppmr vertex correction diagrams due to Nh-loops. 



and the ppnn vertex correction reads as 



^(4 Nh) 
OOab 



^(4 Nh) 
■ jOab 



{k,k,q) 
{k,k,q) 



k 

i^ab - S3aS3b)j^[^Nh{k + q) + IlNh{k - q) - '2IlNh{k) 
{5ab-53a53b)%{TiNh{k-q)-nNh{k + q)) = f%^,^\kXq) 



r 



(4 Nh) . 



ijab {k,k,q) 



-ig'^iSab - S3aS3b)Sij(p-Nh{k - q) + IlNhik + q) 



(4.14) 



Since the result can be expressed through the pion selfenergy H^hi it is quite obvious 
that the Ward-Takahashi identities ( |4.1| ) and ( [4.2| ) are fulfilled for a pion propagator 
dressed by A^/i-loops. 

The vertex corrections which correspond to the A/i part of the pion selfenergy can be 
obtained in the same way. Here we start from the interaction Lagrangians 



= g%/>Ti^^r-lMYp^ + i^pnTpr + l^,Y/>T.^'\ur, (4.15) 



J^^pttNA = -ig—iljf^i!p,p^ ■ Ts(, 



+ h.c. 



(4.16) 



which one gets from gauging £a and Ct,na (eqs. (|3.3|) and ( p.4|) ). As stated in the previous 
section we include a finite width for the A. Here one has to be careful not to violate the 
Ward-Takahashi identity which relates the pAA vertex to the inverse A propagator. This 
is the reason why we have restricted ourselves to a constant A width which drops out if 
we take the difference between two inverse propagators. 

Repeating the entire calculation for the Ah diagrams, the result is completely analo- 
gous to eqs. ( |4.13|) and (|4.14| ) with the index Nh replaced by A/i. For the pTrvr vertex this 
has already been found by Chanfray and Schuck (eq. (3.7) in ref. [TB|). 

The next problem is to include the short-range Migdal interaction. Since the NNNN, 
NNNA and A^A^AA vertices are momentum independent, the p meson does not couple to 
these vertices. Nevertheless many additional diagrams must be calculated. An example for 
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Figure 7: One contribution to the pion selfenergy and the corresponding pTrvr vertex corrections. 



this class of diagrams is shown in fig. 7. It turns out that relations analogous to eqs. (|4.13|) 
and ( 4.14 ) are valid for each pion selfenergy contribution and the corresponding vertex 
corrections. Hence the index Nh in the eqs. ([4.13|) and ( [4.14|) may be omitted. 

So far the tcNN and ttA^A form factor has not been included. Naively one would 
expect that this form factor multiplies the pvrvr vertex corrections by r.^{k)rT^{k + q) and 
the ppniT vertex corrections by r7r(/ci)r^(A;2). However, this violates the Ward-Takahashi 
identities. The reason is that the ttNN (ttNA) form factor implies that the pirNN 
(pnNA) vertex must be modified. We shall do this in a systematic way as proposed in 
refs. Ill, H. 

Formally the form factor defined in eq. ( p.9| ) can be generated by the "propagator" 
l/(— fc^ — A^) of a heavy particle carrying the quantum numbers of a pion (therefore 
called "heavy pion"), which is inserted between the physical pion propagator and the 
ttNN (tcNA) vertex (see fig. 8). Note that we have omitted the /cq in the heavy-pion 
propagator, because depends only on the three momentum k. In order to get the 
correct normalization of the form factor, the "vertex" where the pion is converted into 
the heavy pion must be assigned a factor zA^. 

In this formalism the pion selfenergy S' corresponds to the diagram shown in the 
left part of fig. 9. Here S again denotes the selfenergy without any form factor. Now, 
according to our general prescription, the pvrTr-vertex correction can be constructed by 
coupling the p meson to this diagram in all possible ways. This leads to the three di- 
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Figure 8: Insertion of the "heavy pion propagator" (dotted line) to generate the ttNN and ttNA 
monopole form factor. 





4 y 




Figure 9: Pion selfenergy (left) and pim vertex corrections with wNN (ttNA) form factor 
generated by the "heavy pion". 

agrams which are shown in fig. 9. The first one corresponds to the naive expectation 
that the vertex correction with formf actor is just the vertex correction calculated without 
formf actor multiplied by r^(A;)r^(^ + g). However, there are two more diagrams where 
the p meson directly couples to the "heavy pion" and which are important for preserving 

gauge invariance. 

The vertex of the p meson coupling to the heavy pion is almost identical to the pmr 
vertex. However, since we have omitted the in the heavy pion propagator, we must also 
set the p — component of the vertex equal to zero in order to satisfy the Ward-Takahashi 
identity. Thus the two extra diagrams of fig. 9 contribute to the spatial components of 
the vertex correction only. The final result is 

fS(*.«) = (C>(fc.'i)-«^3.>(2*^+.;).( /^y_^^l^, +^))r.(g)r,(i;+g-) , (4.i7) 

where we introduced the following short-hand notation: for p — the term with the index 
i should be dropped, otherwise i = p. 

In the same way we can construct the ppTrvr-vertex correction with a form factor. 
Again starting from the pion selfenergy as shown in fig. 9 and coupling two p mesons to 
this diagram we find 13 diagrams to evaluate. The result is 
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A2 + (fc - g)2 
(2fc + q)i 



A^ + ik + q}' 



:f IK-/^, -q) + 



A2 + (A; + g)2 
(2A; - g), 



A2 + (A;-g)2 



{2k-qU2k-q),( J:{k - q) ^^_nk) 



A2 + (fc - g)2 VA2 + (fc - g)2 A2 + A;2 
(2A; + g),(2A; + g),^ S(A; + g) ^ 2- ^^^^ 



A2 + (A; + g)2 VA2 + (A; + g)2 A2 + A;2 

S(A:) 



A2 + A;2 



(4.18) 



The meaning of the index i is the same as in eq. ( 4.17| ). Similarly j = v for v = 1,2, 3, 
whereas the corresponding terms vanish for i/ = 0. 



5 The p Meson in Nuclear Matter 

We are now ready to address the central task of our article which is the description of the 
p meson in dense matter. For this we generalize the vacuum selfenergy of the p meson as 
shown in fig. 1 by replacing the vacuum pion propagators and vertices by the in-medium 
propagators and vertex functions as calculated in the previous two sections: 



1 
2 

+ - 



d^k 

(2^ 
d^k 



(2vr) 



^G^{k)V'^ll{k, q)tGAk + q)Kf^{k + q, -q) 
iG^{k)V'}±{k,k,q) . 



(5.1) 



The corresponding diagrams are shown in fig. 10. It is easy to see that S^,^ fulfills the 
transversality condition, eq. (^.51), using the Ward-Takahashi identities, eqs. ( [4.1| ) 
and (|4.2|) one gets 
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c2 ^3ab 



d^k 



{G^k + g)r:S(fc + q, -q) - G^kKfJ^k, -q)) (5.2) 



Since a Pauli-Villars regularization is used we are allowed to substitute k — > k + q in the 
second term and the entire expression vanishes. (Obviously this would not be true if the 
integral was regularized by a form factor.) 



The transversality of S^j, implies its general structure to be 

S^,(g) = ET{q){PT),. + T.L{q){PL)^u , 
with the transverse and longitudinal projection operators 







p = or z/ = 
p,ue {1,2,3} 



(5.3) 



(5.4) 
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Figure 10: p-meson self energy in nuclear matter. 

They are both four-dimensionally transverse, but Pt is three- dimensionally transverse, 
whereas Pl is three- dimensionally longitudinal. Since in matter Lorentz invariance is not 
manifest because there is a preferred frame of reference in which the matter is at rest, the 
two functions and depend on g° and |g| separately. 

Similarly, the in-medium p propagator can be expressed via the projection operators 
Pt and Pl as 



cri?) = GMQ)pr+GpL{q)pr + 



— (mp"'*)^ — Er(g') q'^ — {mf^y — T,L{q) {mp^Yq^ 

We can choose the frame of reference such that q points in the direction of the 2;-axis. 
In this case St = Sn and = (g^/gQ)S33, i.e. it is sufficient to work with the spatial 
components of T^p^. Therefore we have to insert the spatial components of the vertex 
functions calculated in sec. | into eq. (|5.1|). 

First the fco-integration is performed. As in eq. ( |4.3| ) we split the vertex functions 
into a bare part and a vertex correction. The spatial components of the bare vertices do 
not depend on k^. Hence, for these terms the integrands depend on /cq through the pion 
propagators only, i.e. we have to evaluate the expressions 



and 



h = - I dk^ G^{k)G^{k + q) (5.6) 

TT J 

h = - I dk^GJk) . (5.7) 

77 J 

The vertex corrections depend on ko through the pion selfenergies 11'. This leads to 
additional integrals which contain pion propagators and selfenergies, for instance 

h = - I dk^ I\'{k)GJk + q) , (5.8) 

TT J 

which arises from a vertex correction to the ppyrvr- vertex (see fig. 11). Altogether we have 



to evaluate seven fco-i^tegrals /m, ttl = 1,...,7 which are listed in appendix |A.1| . The 
spatial components of S^,^ can then be written as 



= Y E / ft\k, q) Uk, q, go) (5.9) 



m=l 



(2vr) 
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Figure 11: Diagram representing the 1% terms in eq. (5.i). 



with purely real functions fiT'\k,q). The explicit expressions are also given in ap- 



pendix |A.l 



As mentioned earlier, the computation of Sj,- becomes quite involved if we use the exact 



pion selfenergies and propagators as defined in eqs. ( |3.5| ) - ( |3.11| ). Since these functions 
are obtained numerically, also the fco-integrations for the functions Ji to have to be 
performed numerically. Finally they have to be integrated over \k\ and one angle. 

Therefore we begin with the 3-level model introduced in the last paragraph of sec. ^ 
and will come back to the exact case later. In the 3-level model the selfenergies and 



propagators are given analytically (see eqs. (|3.12|) - (|3.16|) ) and also allow for an analytical 



evaluation of the integrals /i to Ij. The results are listed in appendix [A. 2 . The remaining 



integrations in eq. ( |5.9| ), i.e. the integration over the angle 'd between k and q and the 
integration over are performed numerically. 

Results for the p-meson selfenergy as a function of the invariant mass M = \/q^ are 
shown in fig. 12 for different densities g and momenta The real parts of Ex, and 
are displayed in the upper three panels. As discussed in ref. Jl^, in matter, they do no 



longer vanish at M = because of screening effects. The corresponding imaginary parts 
are shown in the lower three panels of fig. 12. Whereas in vacuum Im S^i, vanishes below 
the two-pion threshold, M = 2m^, in matter there are new decay channels, such as Nh 
and Ah states, which allow for finite imaginary parts at lower invariant masses and even 
for M = 0. On the other hand, the imaginary part of the p- meson selfenergy should 
vanish for g° — > 0, independently of q. In our model this is slightly violated because of 
the constant width of the A. An energy- and momentum independent A width means 
that there is a finite probability to find a A with the mass of a nucleon or even below. 
This leads to (unphysical) decay channels of the p meson with zero energy which are 
responsible for the non-vanishing selfenergy at q^ = 0. However, the effect is very small. 
For instance, for Q = 2qq and |g| = g° = we find Im S = lO^^GeV^ which is not visible 
in fig. 12. 

In fig. 13 the longitudinal and transverse parts of the p propagator are shown as a 
function of M and for the same densities and momenta as the selfenergies in fig. 12. As 
expected, the width of the p meson increases with density. Defining the mass of the p 
meson by the zero of the real part in the propagator ( "pole mass" ) , we find an increase with 
density. However, because of the strong broadening, at least q = 2go, the quasiparticle 
concept breaks down. At this density one can also see that the spectral function, i.e. the 
imaginary part of the propagator, develops two maxima, one above and one below the 
vacuum p-meson peak. For q = this was already found in refs. |T3|, At finite q, 
however, the transverse and the longitudinal propagator behave rather differently. 
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q = 0.5 GeV 



q = 0.5 GeV 




Figure 12: From left to right: S = = T,^ for \q\ - 
densities g = (solid line), g = Qo (long dashes) and g 



have been performed using the parameters of refs. ^ 
9i2 = 922 = 



= 0, Et and for \q] = 0.5 GeV at 
= 2^0 (short dashes). The calculations 
m/ (OV = 1200 MeV, = 0.8 and 
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Figure 13: From left to right: Gp = Gpx = Gpi for \q\ - 

densities q = Q (solid line), g = Qo (long dashes) and g 



have been performed using the parameters of refs. ^ 
9i2 = 922 = 



0, GpT and GpL for \q] = 0.5 GeV at 
- 2 go (short dashes). The calculations 
^ (A = 1200 MeV, g[^ = 0.8 and 



18 




Figure 14: Imaginary part of the transverse and longitudinal p propagator Gpx and Gpi at 
g = 2qq as a function of the invariant mass M and the momentum \q\. The calculations 
have been performed using the parameters of refs. j71| , (A = 1200 MeV, g'n = 0.8 and 
9i2 = 922 = O-SJ- 

This can be seen even better in fig. 14 where the imaginary parts of GpT and Gpi at 
g = 2qq are plotted as functions of |g| and M. The two branches we have seen in fig. 13 
are clearly visible in both functions. Note that the lower branch of hxvGpT comes down 
to lower invariant masses with increasing |g| and finally reaches M = 0. This will become 
important later on. 

For the interpretation we take the imaginary part of the propagator: 



Im ST,L(g) 

(g2 - [m';:'f - ReST,L(g))2 + (ImET,L(g))^ 



Im GptA<i) = , . . (0),, ^ ^ L . ^ . ... • (5-10) 



Approximately, the maxima of Im GpT,L correspond to maxima of Im T,t,l or to minima 
of the denominator of eq. ( |5.1CI|) . In vacuum only the latter are present. With increasing 
density the vacuum peaks are pushed to higher as a consequence of an increased real 
part of the selfenergy. This explains the branches of Im GpT and Im GpL at higher 
M. The lower branches correspond to Im Y\t,l and thus, according to eq. ( |5.9|) , to the 
imaginary parts of the functions Ji to I-j. 

From the explicit expressions given in appendix |A.2| (eqs. ( |A.12| ) - (|A.18|) ) and ne- 
glecting the finite A width for the moment we find: 

3 ^ 

Im/„(fc,g,go) = XI V5ir^(^' ^ '^(^o m = 1, 4 , (5.11) 
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Im/6(^,g,go) = E E 4'^^' ^) %o - W + ^.(^ + ^))') (5-12) 
i=i j=i 

and 

— * — * 

Im/5(fc,g,go) = lml7{k,q,qo) = 0, (5.13) 

where (p\^\k,q) are some real, well-behaved functions. Including the A width basically 
leads to some "smoothing" of the (5-functions. 

The structure of the expressions for Im/i to Im/4 suggests the following physical 
interpretation: in matter the decay of the p meson into two pions is replaced by the decay 
into two longitudinal quasiparticles with dispersion relations Ui and ujj (see eq. (|3.16|) ), i.e. 
{Nh)L, TT or {Ah)L quasiparticles. The situation is slightly different for Jg which originates 
from a vertex correction to the ppvrvr-four point vertex (see fig. 11). To this diagram also 
transverse Nh and Ah excitations contribute which have the dispersion relations Qi and 
0,2- In fact, as explained in appendix |A.1| , the decay into two longitudinal modes can 



be separated out and ImJg describes the decay of the p meson into one longitudinal 
quasiparticle and one transverse Nh or Ah excitation. 

Due to their small strength, 5*1, the decay channels involving one or two {Nh)^ quasi- 
particles are relatively unimportant. The decay into two pion quasiparticles is suppressed 
by vertex corrections, as already shown in refs. |13|, Therefore the most important 
contributions to Im/i to ImJ4 come from the decay p Ti{Ah)L. Similarly Im/g is 
dominated by the decay p n{Ah)T- 

In order to identify these two modes in the spectral function we realize that the pionic 
branch has most of its strength at low momenta (see fig. 3). Therefore a p meson with 
momentum q should predominantly decay into a pion with momentum k ^ and a Ah 
pair with momentum close to q. The corresponding invariant masses 



Mp^^(^Ah)r ~ \l {m^ + n2{q)y - , (5.14) 

are shown in fig. 15 for g = 2qq. They are in good qualitative agreement with the 
structures seen in the imaginary part of the p propagator: We can identify the lower 
branch of ImGp^ with the decay mode p —>■ ^{Ah)^ and the lower branch of hxvGpT with 
the decay mode p 7T{Ah)T- Therefore this structure of ImGpT is mainly a consequence 
of the vertex correction in fig. 11. 

All results presented so far have been calculated within the 3-level model, i.e. the 
Fermi motion of the nucleons has been neglected. In the final part of this section we 
want to test the quality of this approximation. Using spectral representations for the 
pion propagator and the selfenergy 11', 



^ , ^lmll'(uj,k) 
n'(k,,k) = (5.16) 
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Figure 15: Invariant mass M of the p meson decaying into a transverse (solid line) or longitu- 
dinal (dashed line) Ah pair, resulting from the approximate formulae eq. (jS.l^ ). 

the imaginary parts of the integrals Ji to J4 and Iq reduce to integrals from to e.g. 

Im/gf/c, g, go) = / dtu lmIl'{ij,k)lmG.,r{qo — uj,k + q) , (5.17) 

TT JO 

and similar expressions for Im Ji to Im I4, while Im J5 and Im vanish. From the structure 
of these integrals one sees immediately that the imaginary part of the p-meson selfenergy 
vanishes for go 0, as it should. As discussed above, because of the constant A width 
this was not exactly true for the results shown before. However, the dispersion relations, 
eqs. ( p.l5|) and ( ^.16| ), are also slightly violated by the constant A width and the two 



effects cancel to some extend giving the expressions like eq. ( 5.17 ) the correct behavior 
for go — > 0. 

The expressions for Im/^ and the remaining angular and momentum integration in 
eq. ( |5.9| ) can be evaluated numerically. Some results for the imaginary part of the p- 
meson selfenergy at g = 2qq are shown in fig. 16. Obviously the difference to the results 
obtained within the three level model is small. This is somewhat surprising since the 
3-level model is known to be a quite poor approximation to the pion selfenergy, at least 
in the domain of the A^/i-excitations. However, after integrating over momenta, details of 
the pion spectral function are "washed out" and only gross features are important which 
are correctly reproduced by the 3-level model. Of course the approximation becomes even 
better at lower densities where the effect of Fermi motion is less important. 



6 Applications: Photoabsorption and Dilepton Pro- 
duction 

In this section we apply our model to the calculation of electromagnetic processes, focusing 
on the effect of the 3-momentum dependence of the pionic part of the p-meson selfenergy. 
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q=0 



q=0.5 GeV 



q=0.5 GeV 




Figure 16: From left to right: Imaginary part of T, = St = El for \q\ = 0, St and St for 
|g| = 0.5 GeV, at density q = 2qq. Solid line: full calculation, dashes: 3-level model. 

We begin with the calculation of the total absorption cross section of real photons on 
nucleons and nuclei p3. In nuclear matter the inclusive cross section per nucleon for a 



photon with momentum q and polarization vector en{q, A) is given by 

^ = l^T.\^Mm\no)\<^o)\\2nrs{pf-q) . (e.i) 

In this expression |/) is the unobserved final hadronic state, |$o) specifies the ground 
state of nuclear matter, and denotes the electromagnetic current operator. 
This expression can be related to the current-current correlation function 

J/^'^(g) = -i j d^xe'^-''{%\T(^.P{x)r{<d))\^Q) . (6.2) 

Inserting a complete set of energy and momentum eigenstates between the two current 
operators and taking the imaginary part, we find 

^ = --^e^(g,A)e.(g,A)ImJ'^^(g) . (6.3) 
A Q\q\ 

In the VDM is a linear combination of the p, uj and field operators. Neglecting the 
isoscalar contributions, the current field identity reads 

jNx) = ^ P ' pf'(x) , (6.4) 
9 

i.e. J^^i^q) is proportional to the p propagator G^''{q). Because of the transverse polar- 
ization of the photon only Gp'^{q) contributes, and we finally obtain 

1 e' {mfY ImST(g) 
A Q\q\ 9^ |(mf )2 + ST(g)|2 ' 
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It is obvious that the 3-momentum dependence of the p-meson selfenergy is necessary for 
the calculation of photoabsorption cross sections within the VDM: Since M = for real 
photons, |g| is the only variable the cross section can depend on. 

For very small densities, eq. ( |6.5|) describes the photoabsorption cross section of a 
single nucleon. Since the p- meson selfenergy in vacuum vanishes for = 0, the limit 
— * is well defined and can be written as 

lim = — 4t ^ lim • (6.6) 

e^o A \q\ e^o Q 

In this limit, interferences between different contributions to S-p vanish and only contri- 
butions to T^T linear in g survive, i.e. diagrams containing exactly one Nh or A/i loop. 
Of course, describing photoabsorption on a single nucleon, the creation of a hole is more 
conveniently interpreted as the destruction of the initial nucleon. 

Since the cross section is proportional to the imaginary part of Sy, the processes which 
contribute correspond to cuts through the transverse part of the p-meson selfenergy. Our 
model, as described in the previous section, leads to the non-resonant background of the 
photoabsorption cross section, e.g. 7 — >■ vr -|- Nh or 7 — >■ vr -|- Ah. The threshold for such 
processes on a single nucleon at rest is |g| > 150 MeV. 

Our results for the non-resonant background of the photoabsorption cross section of 
the nucleon are shown in the left panel of fig. 17. For comparison we also display the 



experimental data, averaged over protons and neutrons The short-dashed line 



corresponds to a form factor cutoff A = 1200 MeV as used in the Bonn potential 
and in refs. |]lT], |12|. In this case, the background alone exceeds the data for go > 800 MeV, 



revealing the cutoff as too large for the present purposes. The solid line was obtained 
with A = 550 MeV, which does not contradict the data. 

Our results are in line with the well-known fact that for many applications a much 
softer form factor is needed than the one used in the Bonn potential (e.g. vrA^ scattering in 
the Chew-Low model Knn ~ 520 MeV, 7r° photoproduction l^: A^^va = 300 MeV, 
deep inelastic lepton scattering Iy-nNN{/\) < 500 MeV). In the cloudy bag model for 



the nucleon ||3^ the form factor has the form 3ji(|A;|i?)/(|fc|-R), which for R = 0.8 fm 
corresponds to a monopole with A^7VAr(A) = VT^/-R ~ 780 MeV This value is similar 
to a result from lattice QCD H), Knn = (750 ± 140) MeV. 

In fig. 17 we also show a decomposition of the background into the processes 7A^ ttN 
(long-dashed line) and 7A^ vrA (dashed-dotted line). Since the A resonance decays into 
7cN (implicitly contained in our model through the constant width of the A) the dashed- 
dotted line can be identified with the two-pion background while the long-dashed line 
corresponds to the single-pion background. The latter is somewhat larger than the single- 



pion background found by Effenberger et al. [40| in an analysis of the 7A^ —>■ irN data. 



This might indicate that our formfactor should be reduced even further . 

^The analysis of pion induced p-meson production (nN pN), which was recently proposed by 
Friman as another test reaction, seems to give a similar result: In order to describe the data we 
would have to reduce the cutoff to about 300 MeV. However, in view of the rather difficult separation of 
the data from competing ttN ttttN processes, further careful studies are required before one can draw 
firm conclusions. See also ref. Es], where the photoabsorption cross section is shown for A = 310 MeV. 
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Figure 17: Left panel: Photoah sorption cross section of the nucleon. The data are the averages of 
the total cross sections for protons and neutrons |^/ . The theoretical calculations correspond 
to the non-resonant background only. They were obtained from eq. \6.^ with g = 0.052 qq and 
cutoff parameters A = 1200 MeV (short-dashed line) and A = 550 MeV (all other lines). The 
background represented by the solid line is also separated into the processes ttN (long- 

dashed line) and vrA (dashed- dotted line). Right panel: Photoab sorption cross section per 

nucleon of uranium. The data are averages of the data from refs. J^l 

550 MeV, g[^ 



the background has 
0.6 and g[2 = 522 = 0.25. 



been calculated using eq. (6J.) with g = 0.8 go, A 
The dashed line was obtained within the 3-level model. The solid line corresponds to an improved 
calculation with the imaginary part of the p-meson selfenergy calculated with Fermi motion while 
for the real part the result of the 3-level model was taken. 
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For a realistic description of the photoabsorption resonance-hole excitations, like 7 
Ah or 7 — > N*{1520)h, have to be taken into account in addition to the non-resonant 
background. The resonant contributions to the p-meson selfenergy have already been 



worked out in a momentum dependent way in ref. |12| and can therefore easily be added 



to our model. This has been done in ref. and we refer to this publication for more 
details. 

In the right panel of fig. 17 we show our results for the photoabsorption on nuclei 



together with the experimental data for uranium |31, 32, 33]. The calculations have been 
done using A = 550 MeV, g[i = 0.6 and g'12 = 922 = 0.25. As an average density in finite 
nuclei we took q = O.Sqq- The dashed line was obtained within the 3-level model. The 
solid line corresponds to an improved calculation with the imaginary part of the p-meson 
selfenergy calculated as outlined in the end of section ^ i.e. including Fermi motion. For 
the real part we kept the result of the 3-level model. Obviously for energies greater than 
~ 150 MeV the differences between the two approximations are small. At lower energies 
the differences in the selfenergy, although still small in absolute units, are enhanced by 
the factor l/\q\ in eq. ( |6.5| ). Therefore the 3-level model should not be trusted in this 
regime. 

In contrast to the photoabsorption on a single nucleon, in nuclear matter there is 
no threshold at |g| ~ 150 MeV. For instance, processes mediated by meson-exchange 
currents, like 7 — tttt — (Nh) (Nh) allow for the absorption of photons with much lower 
energies. Of course our model cannot describe the giant resonances at very low energies, 
but the pionic background contributes considerably to the "dip region" at |g| ~ 150 MeV. 

Again, for a realistic description of the data resonant diagrams have to be taken into 



account. More details about this can be found in ref. p4 . 

Finally we turn to the calculation of dilepton rates in hot nuclear matter. In the VDM 
the rate of e~^e~ pairs with total 4-momentum q produced per volume d^x is given by p 



dN^^^ ^i!^,-ImG« (q-T) (6 7) 

Here a = e^/(47r) is the fine structure constant, M = \fq^ the invariant mass of the 
dilepton pair, and T the temperature of the hadronic matter. G^{q; T) is the retarded p 
propagator at temperature T and baryon density g. Thus eq. ( |6.7| ) depends on temper- 
ature through the Bose factor l/(e'^ ^'^ — 1) and through the retarded propagator. Since 
in sec. ^ the p propagator was calculated for T = only we approximate 

Im Gf{q; T) ^ Im T = 0) , (6.8) 

but keep the temperature dependence through the Bose factor. (Note that at T = the 
imaginary parts of the retarded propagator and the time-ordered one Gp are equal.) 

In fig. 18 the dilepton rates for T = 150 MeV and Q = Qq are shown. In the calculations 
we used the parameters obtained from the fit to the photoabsorption spectra i.e. 
A = 550 MeV, g[^ = 0.6 and 

9i2 — 922 — 0.25. The solid line corresponds to our model 
as described above. For M ~ 0.2 - 0.6 GeV it is in fair agreement with the 'pion -|- one 



loop' result of Steele et al. However, whereas these authors find almost no medium 



effect at and above the free p-meson mass, in our calculation the vacuum peak is strongly 
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Figure 18: Dilepton rates for g = go o-nd T = 150 MeV calculated from eg. H ) with dif- 
ferent approximations to the retarded p propagator G^{q^,q; T): Gp{q^,q\ T = 0) (solid line), 
Gp{M,q = 0; T = 0) (dashed line), and Gp°''^{q^,q; T = 0) (dotted line). The parameters are 
A = 550 MeV, g'^ = 0.6 and g'^^ = 522 = 0.25. 



reduced and we find enhanced rates at higher invariant masses. (For comparison see the 
dotted hne which has been calculated using the p propagator in vacuum.) 

We also show the result of a calculation neglecting the q dependence, i.e. using the 
approximation of refs. [11, 12|, ImG'p(g°, g) ~ Im Gp(M,g = 0) (dashed line). Whereas 
the momentum dependence does not change very much for larger invariant masses, for M 
below ~ 300 MeV the results including the q dependence of Gp (solid line) are significantly 
enhanced as compared to the calculation without q dependence (dashed line). This can 
be explained as follows: As discussed in sec. I, the p meson strongly mixes with 7r(A/;,)T 
states and the invariant mass of these states decreases with increasing q and finally reaches 
zero (see figs. 14 and 15). Because of the factor 1/M^ in eq. ( p.7|) this leads to a strong 
enhancement of the dilepton rate at small M. (Of course, the divergence at M = is 
only an artifact of neglecting the electron mass in the derivation of eq. (|6.7| ).) 

Again, in a realistic calculation resonance contributions of the type p Bh have to 
be taken into account 11121 ITS 



7 Summary and Conclusions 

We have extended hadronic models for p-meson propagation in cold nuclear matter in- 
volving medium modifications of the pion propagator to include finite 3-momenta. The 
starting point was the p propagator in vacuum with a bare propagator dressed by tttt- 
loops. In matter, the pions are renormalized through particle-hole and A-hole excitations 
correlated by Migdal parameters. We use a 3-momentum dependent monopole form factor 
to cut off higher momenta at the ttNN- and vrA^A-vertices. For numerical convenience, in 
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most of our calculations the dressed pion propagator was approximated by the so-called 
"3-level model" which is obtained from the exact propagator by neglecting the Fermi mo- 
tion of the nucleons. This approximation was found to be excellent for our purposes (see 
figs. 16 and 17). 

Because of the medium modifications of its pion cloud, the properties of the p-meson 
become modified as well in matter. Guided by Ward-Takahashi identities, we have calcu- 
lated vertex corrections to the pTrvr- and ppTrvr-vertices in order to preserve gauge invari- 
ance. For the resulting in-medium p propagator these vertex corrections turned out to be 
particularly important to shift part of the strength down to the regime of low invariant 
masses. The main effect could be attributed to the decay of the p into a pion and a 
spin-transverse A-hole excitation. For |^ = this branch was already found by Chanfray 
and Schuck |T^ and by Herrmann et al. [|T^. However, with increasing \q\ it comes down 
to even lower invariant masses and eventually crosses the M = line. 

As one application of the model we calculated the non-resonant background for pho- 
toabsorption on nucleons and nuclei. In order to be consistent with the experimental 
data, the cutoff parameter which enters the ttNN- and vrA^A-form factor had to be re- 
duced from 1200 MeV as used in earlier calculations [O, |12[ to 550 MeV. This example 



shows the importance of constraining the model parameters from experimental data. As 
discussed in section ^ there are indications from other observables, like more exclusive 
photoabsorption data (see ref. [0) or pion-induced p-meson production |^ that a further 
reduction of the cutoff might be necessary. This is presently under investigation. 

Finally we have studied the consequences of the 3-momentum dependence of the p- 
meson propagator for dilepton production rates in hot and dense matter. Compared to 
a calculation where the |g| -dependence was neglected we found almost no difference for 
invariant masses above ~ 400 MeV but a strong enhancement of the rates at lower M. 
This is related to the 7r(A/;,)T-branch as discussed above. 

In the present article we have concentrated on the pion-loop contribution to the p- 
meson selfenergy. However, realistic calculations require to account for contributions 
arising from direct resonance formation p —>■ Bh [jl2|, |T6|. In fact, the medium modifi- 
cations of the p propagator associated with the pionic part of the selfenergy are substan- 
tially decreased by the use of much softer irNN and ttA^A form factors. This increases the 
relative importance of the resonances as compared to earlier versions of the model where 
a larger cutoff was used |jl2[ . As long as interferences induced by explicit resonance decays 



are neglected, the resonance-hole diagrams can be evaluated separately and simply added 
to the pionic selfenergy contributions |2i|. Ultimately, one should take into account the 
resonance decays into vrA^ and pA^ in a selfconsistent fashion [l^, which also implies a 
better description of the A width. 
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A p Meson Selfenergy at Finite Density 



A.l General Expressions for the Spatial Components 

As described in sec. ^ the spatial components of the p-meson selfenergy are calculated by 
inserting the spatial components of the vertex functions calculated in sec. I into eq. 
If we split the vertex functions into bare vertices and vertex corrections, as in eq. ( [4.3| ), 
the selfenergy becomes a sum of seven terms: 

-2 " d^'k 



with 



+4 



{ki{K^-P)-qiP){2k + q), 



A^ + {k + q)^ 

^ ^ ihjA^ - P) - q.P)ik,iA^ -ik + q)') + q,A') ^ 
(A2 + (A?+g)2)(A2 + p) ^ 

+2/5 

^ JA2 + P)6^, - {2k + g),(2fc,A2 - q^P) ^ 

(A2 + (fc + g-)2)2 ' 
_4f + + g).(A;,(A2 - P) - q^P] 



-A2 + A;2 (A2 + (A; + g)2)(A2 + p) 
+ (A.l) 



(A.2) 
(A.3) 
(A.4) 
(A.5) 
(A.6) 
(A.7) 
(A.8) 



h = 


nJ 


f dk'^GMGAk + q) , 


h = 


^ J 


f dk^U'{k)G^{k)G^{k + q) , 


h = 




f dk^U'\k)G^{k)G^{k + q) , 


h = 




f dk^ U'{k)G^{k)U'{k + q)G^{k + q) 


h = 




f dk^G^ik) , 


h = 


nJ 


f dk^U'{k)G^{k + q) , 


h = 




f dk°U'{k)G^{k) . 



It is useful to rearrange these terms slightly [|T^: in the integrals I2, I3, I4, and Jy the 



vertex corrections 11' are always coupled to a pion propagator of the same momentum. 
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Therefore, only the spin-longitudinal Nh or Ah excitations contained in 11' can contribute 
to these integrals. This is different for Iq. Here longitudinal and transverse Nh and Ah 
excitations of the vertex correction can contribute, because they are not coupled to a 
pion. However, we can separate longitudinal and transverse Nh and Ah excitations in 
the Iq term, and combine the longitudinal part with the Is term. To that end, we define 
the longitudinal spin-isospin response function, 

U'Lik) = U'{k) + U'{k)k^G^{k)U'{k) = {k^ - ml)U'{k)G^{k) , (A.9) 

and an integral similar to Iq, 

h = - I dk' U'^ik)G^{k + q)^ k% + h ■ (A.IO) 
With this abbreviation we can write 

f (f'k 



^iM = ^Jj^{i^k + qU2k + q)jh 

-P)-q,P){:. 
K^+ik + q)'^ 



(27r)= 

, {k,{K^-P)-q,P){2k + q)^ ^ 

+4 — — — J2 



^ , (fc.(A^ - P) - g.P)(A;,(A^ - k^) - q,P) ^ 

^ ^ (fc,(A2 - P) - q,P){k,{A' - jk + q)') + q,A') ^ 
{A^ + {k + q)^){A^ + P) 

(A2 + (A; + g)2)2 V ' ^2 ) 
^( k%j ^ (2fc + g).(fc,(A2-fc2)-g,fc2) ^ X 

^A2 + ^2 (A2 + (^+g)2)(A2 + fc2) J V 

+ (A.ll) 

A. 2 3-Level Model 

Within the 3-level model the integrals Ii to can be evaluated analytically: 



3 3 



Si{k) Sj{k + q) Ui{k) + Uj{k + q) 



i=i j=i uji{k) ujj{k + q) ql- {uJi{k) + Uj{k + q)Y 
^ j^j^ Qjk) Sjjk + q) uj,{k) + UJ,{k + q) ^^^^^ 

i=i j=i iOi{k) ujj{k + q) ql - {uJi{k) + LOj{k + q)Y 
^ j^j^ {u;f{k)-ool{k))Q{k) S,{k + q) u;,{k)+u;,{k + q} 



uJi{k) ujj{k + q) qi - {oJi{k) + ujj{k + g))2 
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= ^^ ^'.(fc) C,{k + q) Ui{k)+UJj{k + q) ^^^^^ 
i=ij=i uJi{k) ujj{k + q) ql- {uJi{k) + ujj{k + q)Y ' 

/5 = (A.16) 

i=i (^i{k) 

jg = ^^ a^(fc) + + ^^^^^ 

Ir = (A.18) 

i=i uJi{k) 

The abbreviations Ci{k) are defined similarly to the strengths Si{k) (see eq ( |3.16| )) by the 
equation 

(U'G.m ^ + + . (A.19) 

k',-ul{k) kl-uol{k) kl-ul{k) 
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